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Abstract 
Recently released Planck data favor a lower value of the Hubble constant and a higher value of 
the fraction matter density in the standard ACDM model, which are discrepant with some of the 
low-redshift measurements. Within the context of this cosmology, we examine the consistency of 
the estimated values for the Hubble constant and fraction matter density with redshift tomography. 
Using the SNe Ia, Hubble parameter, BAO and reduced CMB data, which are divided into three 


bins, we find no statistical evidence for any tension in the three redshift bins. 
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I. INTRODUCTION 


More than one decade ago it was found that our universe is in an accelerating expansion 
based on the distance measurement of type Ia supernovae (SNe Ia) Al Lh. This observation 
is consistent with other astronomical observations such as Hubble parameter, large scale 
structure and cosmic microwave background radiation (CMB), etc. To explain this acceler- 
ating expansion, one has to introduce the so-called dark energy with negative pressure in the 
general relativity framework, or to modify the general relativity at cosmic scales. Although 
suffered from some theoretical issues, the cosmological constant introduced by Einstein 
himself in 1917 is the most simple and economical candidate for the dark energy. Indeed 
the standard ACDM model turns out to be consistent with several precise astronomical 
observations, such as SNe la l, Wilkinson Microwave Anisotropy Probe (WMAP) mea- 
surements of CMB |6], and baryon acoustic oscillation, etc. If the standard ACDM model 
properly describes our universe, the current Hubble constant Ho and fraction matter density 
Omo Should be consistent with those estimated by different observations made at different 
redshifts. 

However, the recently released Planck data td favor a higher value of Qmo = 0.315+0.017 
and a lower value of Hp = (67.3+1.2) km s~' Mpc™ in the standard six-parameter ACDM 
cosmology, obtained by using Planck+WP, where WP stands for WMAP polarization data. 
These values are in tension with the magnitude-redshift relation for SNe Ia and recent 
direct measurements of Ho, such as the Hubble_space telescope observations of Cepheid 
variables with Hy = (73.8 + 2.4) km s~! Mpc™! [8] and Ho = [74.3 + 1.5(stat.) + 2.1(sys.)] 
km s-' Mpc7! obtained by using a mid-infrared calibration of the Cepheid distance scale 
based on observations at 3.6 um with the Spitzer Space Telescope |9]. Of course, if relax the 
restriction of the standard six-parameter ACDM model, for example, consider the dynamical 
dark energy model or include the dark radiation l, the tension might be alleviated. In 
12|, Hu et al. found that there is another way to alleviate this tension in modified gravity 
models. Furthermore it was reported in that this tension may also be alleviated by if 
one first calibrates the light-curve fitting parameters in the distance estimation in SNe Ia 
observations with the angular diameter distance data of the galaxy clusters, with the help 
of the distance-duality relation. Very recently, Efstathiou reanalyzed the Cepheid data 
and found Hp = (70.6 + 3.3) km s7} Mpc™! based on the NGC 4258 maser distance and 
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Ho = (72.5 + 2.5) km s~' Mpc7! with three distance anchors combined, which alleviates 
the tension compared to the result obtained by Riess et al. g , but the latter still differs 
by 1.90 from the Planck value. In addition, by comparing the eight ultra low redshift SNe 
Ia data (z = 0.0043 to 0.0072) with low redshift data (z < 0.04) from the Union2.1 
compilation and Planck data d. Zhang and Ma found that the present expansion of 
the universe estimated from the low redshift measurements is higher than the one estimated 
from high redshift observations in the ACDM model . In other words, higher redshift 
measurements give a lower value of h, the reduced Hubble constant. 

These discrepancies seemingly imply that the standard ACDM model cannot well describe 
the properties of the universe at all redshift if the major sources of systematic errors of these 
observations have been controlled. In this paper we detect these discrepancies in the ACDM 
model with redshift tomography. We divide the redshift range under consideration into three 
bins and use observation data in each bin to separately constrain the Hubble constant and 
fraction matter density in the ACDM model. In the literature the redshift tomography is 
often used to see the dynamical property of dark energy by piecewise parametrization of the 
equation of state of the dark energy. Here our goal is to see the consistency of the ACDM 
model at different redshifts, therefore we focus on the ACDM model. The data sets we use 
here include the Union2.1 SNe Ia data [15], 19 Hubble parameter H(z) data , Baryon 
Acoustic Oscillation (BAO) data measured by the 6 degree Field Galaxy Survey (6dFGS), 
SDSS DR7, SDSS DR9 and WiggleZ surveys, reduced nine-year WMAP data (WMAP9) 
and reduced Planck data both based on the flat ACDM model. 

The paper is organized as follows. In section [LI] we describe the redshift tomography 
method and observational data. In section we show the results of different combination 
of data sets to constrain the base ACDM model based on the SNe Ia data and the redshift 


tomography analysis. The results are summarized in section 


Il. METHOD AND DATA 


In a spatially flat Friedmann-Robertson-Walker universe, the Hubble parameter is given 


by the Friedmann equation 


H?(z) = HG [Qro(1 + 2)* + Qamo(1 + 2)? + Qo(1 + 2)? + (1 — Rmo — Qrg)] , (1) 
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for the ACDM model, where the redshift z is defined by (1 + z) = 1/a, and Qro, Qumo 
and Qe are the present values of the fraction energy density for radiation, dark matter and 
baryon matter, respectively. The latter two are often written as the total matter lh 


Omo = Qo +Qamo. The radiation density is the sum of photons and relativistic neutrinos 
Qro = QP (1 + 0.2271.Ne fs), (2) 


where Nest = 3.046 is the effective number of neutrino species in the Standard Model Lod, 
and QP = 2.469 x 10-5h-? for Toms = 2.725K (h = Ho/100 km s~! Mpc7?). 

We focus on constraints on the Hubble constant and the fraction matter density in the 
context of the ACDM cosmology from the low-redshift observational data including Union2.1 
SNe Ia sample, Hubble parameter and BAO data, in combination with the high-redshift 
CMB measurements. We adopt a redshift tomography method to examine the flat ACDM 
model. Since the SNe Ia data cannot alone constrain the ACDM model very well, it could 
be even worse in each redshift bin because of the decreasing of data points (so do the Hubble 
parameters), we then divide the redshift into three bins so that the BAO data can distribute 
uniformly in the first two bins, while the CMB data are in the third bin. As a result, these 
data are divided into three combinations in the following redshift bins: 0 — 0.28, 0.28 — 0.73 
and > 0.73. The distribution of data is listed in Table [1] To see the difference, we will use 
WMAP%9 and Planck data separately. 


EF "h 5äF, DRva, REAO 


Table 1: Distribution of SNe Ia, Hubble, BAO, CMB data in three redshift bins. 


The best-fit values of Qmo and h, and their 68% and 95% confidence level (CL) errors are 
obtained by performing the Markov Chain Monte Carlo analysis in the multidimensional 
parameter space in a Bayesian framework. Since the Hubble constant is completely degen- 
erate with the absolute magnitude of SNe Ia, SNe Ia data are not sensitive to the Hubble 
constant. Therefore, in our analysis we marginalize analytically over the Hubble constant 
when the SNe Ia data are concerned. Moreover, note that the fraction baryon energy density 


Qvo is involved in the likelihood for the BAO and CMB data. 
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A. Type Ia Supernovae 


The SNe Ia data set is an important tool to understand the evolution of the universe. 
In this work, we adopt the Union2.1 compilation , containing 580 SNe Ia data over the 
redshift range 0.015 < z < 1.414. The chisquare is defined as as 

N oos 
[u (zi) — (2)? 


2 
FA = n (3) 
= i=l Tay) 


where N is the data number in the redshift interval we are interested in, 1°°%(z) is the 
measured distance modulus from the data and p"(z) is the theoretical distance modulus, 
defined as 

u (z) = 5logio dr + Ho, Ho = 42.384 — 5 logio h. (4) 


The luminosity distance is 
dz 


dilz) = (1+2) [ EG)’ (5) 


where E(z) = H(z)/Hpo. We can eliminate the nuisance parameter jig by expanding x? with 


respect to uo [21 : 


Xgn = A+ 2B uo + Cys, (6) 
where 

N 

re 5 [u (zi; po = 0) — u% (2)? 

= a 

=i ogy (zi) 
AL pth (za; po = 0) — u* (z) 

B = ` — 2 ? (7) 
ri ogy (zi) 


The 2, has a minimum as 


Kew = A-B/C , (8) 


which is independent of jig. This technique is equivalent to performing a uniform marginal- 


ization over Ho . We will adopt %Zy as the goodness of fitting instead of y2)y. 
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B. Observational Hubble parameter (HUB) 


The observational Hubble parameter can be obtained by using the differential ages of 


passively evolving galaxies as 


1 Az 
H == 9 
1+z At (9) 


We use 19 observational Hubble data over the redshift range 0.07 < z < 2.3, which con- 
tain 11 observational Hubble data obtained from the differential ages of passively evolving 
galaxies , and 8 H(z) data at eight different redshifts obtained from the differential 
spectroscropic evolution of early type galaxies as a function of redshift |19]. The chisqure is 
defined as : 


ue _ 2 Palo Aei , (10) 


where H;;(z) and Hos(z) are the theoretical and observed values of Hubble parameter, and 


oy denotes the error bar of observed data. 


C. Baryon Acoustic Oscillation 


As aruler to measure the distance-redshift relation, Baryon Acoustic Oscillation provides 
an efficient method for measuring the expansion history of the universe by using features in 
the cluster of galaxies with large scale surveys. Here we use the results from the following 
five BAO surveys: the 6dF Galaxy Survey, SDSS DR7, SDSS DR9, WiggleZ measurements 


and the radial BAO measurement. 


1. 6dF Galary Survey 


The 6dFGS BAO detection allows us to constrain the distance-redshift relation at Zeff = 
0.106 22), The low effective redshift of 6dFGS makes it a competitive and independent 
alternative to Cepheids and low redshift supernovae in constraining the Hubble constant. 
They achieved a measurement of the distance ratio 


Ts(za) 


— Aa O = 0.336 + 0.015 ti 
Dy(z = 0.106) i P 
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where r,(za) is the comoving sound horizon at the baryon drag epoch when baryons became 
dynamically decoupled from photons. The redshift zg is well approximated by 


1291 (Qmoh?)0251 


EE P O [1 + bi (Qh) *], (12) 


Žd 


where 


bi = 0.313(Qmok OS [1 + 0.607 (Qmo)? 6], a3) 
by = 0.238 (moh). 
The effective “volume” distance Dy is a combination of the angular-diameter distance D 4(z) 
and the Hubble parameter H(z), 
= dr z 1 
Dy(z) = [ al 
o A(x)’ H(z) 


(14) 
=| 2) Dale) ay). 
The xZqp is given by ; 
Xar — Kalaa) Py 0 = 0.336] ; (15) 


2. SDSS DR7 


The joint analysis of the 2-degree Field Galaxy Redshift Survey data and the Sloan Digital 
Sky Survey Data Release 7 data gives the distance ratio at z = 0.2 and z = 0.35 [24]: 


_ Ts(24) L 0.1905 + 0.0061, 
Dy(z = 0.2) (16) 
__ "s(2a) _ 9.1997 + 0.0036. 
Dy(z = 0.35) 
When the two data points are in the same redshift bin, we adopt the y7p7 given by 
XDRT = X"V“X, (17) 
where 
rs(Za) 
rs(Za) 
PAA — 0.1097 


and the inverse covariance matrix is 


30124.1 —17226.9 
—17226.9 86976.6 
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On the other hand, when the two data points are in the different redshift bin, their X5 p7 


are respectively given by 


[( 42), — 0.1905]? 


eee 0.0061? i (20) 
; (ets en — 0.1097)? 
ADENS 0.00362 


3. SDSS DR7 reanalysis 


By applying the reconstruction technique to the clustering of galaxies from the SDSS 
DR7 Luminous Red Galaxies sample, and sharpening the BAO feature, Padmanabhan et 
al. obtained the distance ratio at z = 0.35 [26] : 


Ts(za) 


—— 5 _ _ = 0.1126 + 0.0022. 21 
Dy(z = 0.35) 21) 


The Xbpry-re used in the Markov Chain Monte Carlo analysis is 


2 (ees Jin — 0.1126)? 
Xbre = 9.00097 (22) 


Since the SDSS DR7 and SDSS DR7 reanalysis results are based on the same survey and the 
latter gives a higher precision than the former, we include the SDSS DR7 reanalysis data 
when we do the whole redshift analysis but not both together. On the other hand, when we 
do redshift tomography, we may refer to part of the SDSS DR7 data at z = 0.2 and when 
the redshift bin contains z = 0.35, we will use the SDSS DR7 reanalysis data. 


4. SDSS DRI 


The SDSS DR9 measurement at z = 0.57 analyzed by Anderson et al. gives 


rs(2a) 


-Sa O = 9.0732 + 0.0012 23 
Dy (z = 0.57) ' 2) 


which is the most precise determination of the acoustic oscillation scale to date. The 


chisquare is defined as 


(pagan — 0.0732)? 


2 
XDR9 = 0.00122 
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5. The WiggleZ measurements 


The WiggleZ team encodes some shape information on the power spectrum to measure 
the acoustic parameter [28]: 

A(z) = ——————_.. (25) 

The measurements of the baryon acoustic peak at redshifts z = 0.44, 0.6 and 0.73 in the 

galaxy correlation function of the final dataset of the WiggleZ Dark Energy Survey give the 


acoustic parameter: 
A(z = 0.44) = 0.474 + 0.034, 
A(z = 0.60) = 0.442 + 0.020, (26) 
A(z = 0.73) = 0.424 + 0.021. 


The chisquare is defined as 
Xwig = X'V*X, (27) 
where 
A(z = 0.44) sn, — 0.474 
X = | A(z = 0.60), — 0.442] , (28) 
A(z = 0.73)ın — 0.424 
and its inverse covariance matrix is 
1040.3 —807.5 336.8 
V-* = |-807.5 3720.3 —1551.9| . (29) 
336.8 —1551.9 2914.9 


6. Radial BAO 


The radial (line-of-sight) baryon acoustic scale can also be measured by using the SDSS 
data. It is independent from the BAO measurements described above, which are averaged 
over all directions or in the transverse directions. The measured quantity is 

A, (2) = H (z)rs(z4); (30) 
whose values are given by 
A, (0.24) = 0.0407 + 0.0011, 


(31) 
A. (0.43) = 0.0442 + 0.0015. 
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D. Cosmic Microwave Background 


In the CMB measurement, the distance to the last scattering surface can be accurately 
determined from the locations of peaks and troughs of acoustic oscillations. There are two 
quantities: one is the “acoustic scale” 


TDa(2Zx) 


l4=(1+ ee 


(32) 


and the other is the “shift parameter” 


R=4/ 0, tei. (33) 


These quantities can be used to constrain cosmological parameters without need to use the 


full data of WMAP9 ld. Here z, is the redshift at the last scattering surface 
2, = 1048/1 + 0.00124 (Qh?) 73) [1 + g1(Qmoh?) I, (34) 


where 


0.0783 (Oa h 

1 £39.5(Qyoh?)2-763" 
7 0.560 

14 21.1(Qe9h?)2 8" 


gı 
(35) 


g2 


Wang and Wang have obtained the mean values and covariance matrix of 
{R, la, Qnoh?, ns} from WMAP9 and Planck data respectively, based on the ACDM model 
without assuming a flat universe. On the other hand, Shafer and Huterer derived the 
related results about {R, l4, z.} from WMAP9 and Planck data respectively, based on the 
flat wCDM model. For our propose, following and Bi we first extract the mean values 
and covariance matrix of {R, l4, z4} from WMAP9 and Planck data respectively based on a 


flat ACDM model. 


1. WMAP9 


By using the WMAP9 data, we obtain the mean values for {R, l4, 24} as 


(la) = 301.95, (R) = 1.7257, (z,) = 1088.96. (36) 


10 
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Their inverse covariance matrix is 


3.087 15.160  —1.456 
Ci mapo = | 15460 198053: —217.021 
1.456 —217.021 5.552 


The chisquare for the reduced WMAP9 data is defined by 


2 = Try-1 
Xwaapo = X Cyapo*: 


where 
la — 301.95 
X= | R—1.7257 
Zą — 1088.96. 
2. Planck 


By using the Planck data, we obtain the mean values for {R, la, z.} as 
(la) = 301.65, <A) = 1.7500, <2.) = 1090.33. 


Their inverse covariance matrix is 


40.909 —405.455 —0.5443 
Colanck = | —405.455 55662.83 —751.123 
—0.5443 —751.123 14.6187 


The chisquare for the reduced Planck data is defined as 


2 — Try-1 
XPlanck = x C PianckX> 


where 
la — 301.65 
X = | R— 1.7500 
zą — 1090.33. 


II. RESULTS 


(38) 


(40) 


(41) 


(42) 


(43) 


Using the Union2.1 sample in combination with other measurements described in the 


previous section, we give the constraints on the base ACDM model. The best-fit values of 


i 


201709.00161v1 


chinaXiv 


Qmo and h with 68% CL errors are summarized in Table 2] and their likelihoods are shown 
in Figure 

From Table 2] we can see that the SNe Ia data alone favor a lower value of Qmo than the 
SNe Ia data in combination with other datasets. Including the BAO and WMAP9/Planck 
data improves significantly the constraint on Qmo. Including the reduced Planck data gives 
the highest Qmo and the lowest h. We find these estimates of Qmo are consistent with each 
other within lo CL, but are in tension with the results derived by Planck [7]. The estimates 
of h from the HUB, BAO and WMAP9 are compatible with those from Planck, but are 
discrepant with those from fitting the calibrated SNe magnitude-redshift relation n 


data Ow Mh 


SN 0277677 faa | — 


+0.0231 +0.0234 
SN + HUB 0.2834 00227 |0. 7055 10.0248 


DADE RA a a E 
SN + HUB + BAO + W M AP9|0.2877+ 00116 |0-702219 age 
SN + HUB + BAO + Planck 009TA hice 


Table 2: Constraints with lo errors on Qmo and h for the base ACDM cosmology from SNe Ia data 


0.703270 vios 


SN + BAO + WMAP9Ì 


in combination with HUB, BAO, WMAP9 and Planck. 


Using the SN+HUB+BAO+WMAP9/Planck data distributed in three different redshift 
bins, we present the constraints on Qmo and h for the ACDM model in Table The 
corresponding marginalized posterior distributions are shown in Figure 

Our analysis shows that low-redshift observations give a higher value of Qmo, while high- 


redshift observations give a lower one by using the SN+HUB+BAO+WMAP%9 data. How- 


ever, the high redshift z > 0.73 observations with Planck data favor a relatively higher value 


of Qmo, which is inconsistent with the high-redshift value from WMAP9 at about 1.lo CL. 


12 
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———— __ Snla+Hubble 

— __ Snla+Hubble+BAO 
———— __ Snla+Hubble+BAO+WMAP 
0.87 


= Snla+BAO+WMAP 


Snla+BAO+Hubble+PLANCK 


0.6F 


P(h) 


0.47 


0.2 


“0.62 0.64 0.66 0.68 0.70 0.72 0.74 0.76 


=  Snla 
\ ——— __ Snla+Hubble 

\ ——— __ Snla+Hubble+BAO 
——._ Sila BAO+ WMAP9 
0.8} 


——= Snla+Hubble+BAO+WMAP9 -4 


me Snla+Hubble+BAO+PLANCK 


0.6F 


P(Qmo) 


0.4F 


0.2 + 


0.0 


0.20 0.25 0.30 0.35 


Fig. 1: Marginalized posterior distributions for Qmo (upper) and h (bottom) from SNe Ia data in 
combination with HUB, BAO WMAP9 and Planck. 


In addition, there are large uncertainties in the estimation of Qmo from the data in the 
redshift range 0 < z < 0.28. From Table B] we find that the data in the mid-redshift range 
0.28 < z < 0.73 favor a lower Hubble constant with a little large uncertainty than the data 
at low and high redshifts. Figure [8] shows the best-fit values of Qmo and h with lo errors 
for the data in three different redshift bins. 


13 
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redshift range h 

0 — 0.28 0.6870) 05 
0.28 — 0.73 0.6677 5.0469 
> 0.73(W MAP) 0:7123 | eae 
> 0.73(Planck) 0.6831 16:0181 
whole(W M AP) 0.287740.0114]0.7022+00110 
whole(Planck) 0.2969+9:0103 |0.6974+9:0082 


Table 3: Constraints with lo errors on Qmo and h for the base ACDM cosmology in three redshift 


bins from the SN+HUB+BAO+WMAP9 and SN+HUB+BAO+Planck. For a comparison, we 


also list the constraints in the whole redshift range. 


In our analysis the Hubble constant is marginalized as a nuisance parameter in the SNe Ia 
likelihood function. Therefore, the constraints on h mainly come from the HUB, BAO and 
WMAP9/Planck data. In Ref. [16] a higher value of Hubble constant is recently obtained 
from measurements of nearby SNe Ia with help of measurements of Cepheid variables, than 
that obtained by Planck. However, our estimates of h from the data in the redshift ranges 
of z < 0.28 and 0.28 < z < 0.73 are lower than the result obtained in [16]. Moreover, the 
high-redshift data (z > 0.73) including the WMAP9 data favor a higher value of h than the 
data in the first two redshift bins. 


IV. CONCLUSIONS 


The estimates of Qmo and h in the base ACDM model should be consistent with each 
other from measurements made in different redshift intervals, if the simplest ACDM model 
completely describes the evolution of our universe and the unknown sources of systematic 
errors of these measurements can be negligible. The recent Planck observations of the 
CMB lead to a Hubble constant of h = 0.673 + 0.012 and a matter density parameter of 
Omo = 0.315 + 0.017 [7], which, however, are different from the low-z measurements. In this 
work, we have studied the consistency of the estimated values for the Hubble constant and 
matter density parameter from different redshift data. 


We have first obtained reduced CMB data for {R, la, 2+} from WMAP9 and Planck data, 


14 


201709.00161v1 


chinaXiv 


SE binl 

— bin2 

— bin3 (WMAP) 
aaa bin3 (PLANCK) 4 
— whole (WMAP) 
whole (PLANCK) 


P(h) 


-- binl 
— bin2 
= bin3 (WMAP) 
=-  bin3 (PLANCK) = 
— whole (WMAP) 


— whole (PLANCK) 


P(Omo) 


Ono 


Fig. 2: Marginalized posterior distributions of Qmo (upper) and h (bottom) for three redshift bins 
of the SN+HUB+BAO+WMAP9/Planck data. 


based on a flat ACDM model. We then have placed constraints on the base ACDM model 
using astrophysical measurements of SNe Ia, Hubble and BAO, in combination with the 
reduced WMAP9/Planck CMB data. We have found that the SNe Ia data alone favor a 
lower value of Qmo and adding the HUB, BAO, and the reduced WMAP9/reduced Planck 
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0.80 — 
0.75 + J 
= 0.70- i | 


0.65 + 7 


0.60 f 1 1 1 ii 1 1 1 1 1 1 1 1 1 1 1 1 1 f 1 1 1 1 1 1 1 1 1 1 1 
0.0 0.5 1.0 1.5 2.0 2.5 3.0 


0.4L | 


Omo 


0.3 i 


0.2L | 


0.0 0.5 1.0 1.5 2.0 2.5 3.0 


Fig. 3: Best-fit values of Qmo (upper) and h (bottom) with 1ø errors in three redshift bins and 
in the whole redshift range. The third and forth error bars in black represent the case containing 


WMAP9 data, while the red ones represent the case containing Planck data. 


data can give a higher one, but it is still in tension with the result reported by Planck. 
Moreover, the estimates of h from the HUB, BAO and WMAP9 are compatible with those 
from Planck, but_are discrepant with those from fitting the calibrated SNe magnitude- 


redshift relation |8]. There is no any tension on h among three redshift bins, as shown in 


16 
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Fig. 

We have also implemented the redshift tomography analysis in the context of the ACDM 
cosmology with the SNe Ia, HUB, BAO and CMB data. We have found that low-redshift 
observations (z < 0.28) give a higher value of Qmo, as estimated by Planck, while high- 
redshift observations (z > 0.73) with the WMAP9 data give a lower one, which is inconsistent 
with that from the SN+HUB+BAO data in the high-redshift range in combination with the 
Planck data at about 1.lo CL. In addition, the data in the mid-redshift range 0.28 < z < 
0.73 favor a lower Hubble constant. The current data cannot provide statistically significant 


evidence for any tension among the different redshift bins. 
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